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ABSTRACT 

A broadband study of the high redshift blazar S5 0836+7 \{z = 2.172) is presented. 
Multi-frequeney light eurves show multiple episodes of X-ray and 7 -ray flares, while 
optieal-UV fluxes show little variations. During the GeV outburst, the highest 7 -ray 
flux measured is (5.22 ± 1.10) x 10“® ph cm“^s“^ in the range of 0.1—300 GeV, 
whieh eorresponds to an isotropie 7 -ray luminosity of (1.62 ± 0.44) x 1050 gj.g g 
thereby making this as one of the most luminous 7 -ray flare ever observed from any 
blazar. A fast 7 -ray flux rising time of ~3 hours is also notieed whieh is probably 
the first measurement of hour seale variability deteeted from a high redshift {z > 2) 
blazar. The various aetivity states of S5 0836+71 are reprodueed under the assumption 
of single zone leptonic emission model. In all the states, the emission region is loeated 
inside the broad line region, and the optieal-UV radiation is dominated by the aeeretion 
disk emission. The modeling parameters suggests the enhaneement in bulk Lorentz 
faetor as a primary eause of the 7 -ray flare. The high X-ray aetivity with less variable 
7 -ray eounterpart can be due to emission region to be located relatively closer to the 
black hole where the dominating energy density of the disk emission results in higher 
X-ray flux due to inverse-Compton scattering of disk photons. 

Subject headings: galaxies: active — gamma rays: galaxies — quasars: individual (S5 
0836+71) — galaxies: jets 
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1. Introduction 

Blazars are luminous active galactic nuclei (AGN) with powerful relativistic jets aligned close 
to the line of sight to the observer (jet viewing angle 9^ < 1 /T, where F is bulk Lorentz factor). 
Due to their peculiar orientation, the radiation observed from blazars is dominated by the emission 
from relativistic jets (Blandford & Rees 1978) which transport momentum and energy to large 
scales. Blazars include the luminous flat spectrum radio quasars (FSRQs) characterized by strong 
and broad optical emission lines (rest frame equivalent width > 5 A) and BL Lac objects where 
optical lines are weak or absent (Stocke et al. 1991; Stickel et al. 1991). A new classification 
based on the broad line region (BLR) luminosity (Lblr) in units of Eddington luminosity (LEdd) 
is proposed by Ghiselliniet al. (2011). In this scheme, FSRQs are defined as sources with 
-^BLR/-^Edd > 5 X 10“^. In addition to that, both classes share many common properties, such as 
rapid flux and polarization variations (Wagner & Witzel 1995; Andruchow et al. 2005), flat radio 
spectra (a^ < 0.5; Sy oc z/““) at GHz frequencies, and exhibit superluminal patterns at radio 
wavelengths (Jorstad et al. 2005). 

In the widely accepted scenario of blazars, a single population of high energy electrons 
present in the jet, emits over the entire electromagnetic spectrum via synchrotron and inverse 
Compton (IC) mechanisms, the former being relevant at lower energies and the latter dominating 
at higher energies. Consequently, the spectral energy distribution (SED) of blazars exhibit a 
double hump shape, with low energy synchrotron peak lying in between infrared to X-rays and 
the high energy IC hump extending up to GeV/TeV 7 -rays. Possible sources of seed photons for 
IC scattering could be internal (synchrotron self Compton or SSC, Konigl 1981; Marscher & Gear 
1985; Ghisellini & Maraschi 1989) or external (external Compton or EC, Begelman & Sikora 
1987; Melia & Konigl 1989; Dermer et al. 1992) to the jet. The plausible reservoir of seed 
photons for EC can be the accretion disk (Dermer & Schlickeiser 1993; Bottcher et al. 1997), the 
BER (Sikora et al. 1994; Ghisellini & Madau 1996), and the dusty torus (Blazejowski et al. 2000). 
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Moreover, the presence of the high energy peak is also attributed to hadronic processes initiated 
by relativistic protons co-accelerated with the electrons (e.g., Miicke et al. 2003; Bottcher et al. 
2013). Besides non-thermal jet radiation, the thermal emission from the accretion disk is also 
observed in many blazars (see e.g., Ghisellini et al. 2010a). Blazars are found to follow a trend, 
so-called ‘blazar sequence’, with more luminous sources (generally FSRQs) being more 7 -ray 
dominated, and have both synchrotron and IC peaks located at lower energies than their fainter 
(mostly BL Lac objects) analogues (Fossati et al. 1998). However, Giommi et al. (2012) have 
recently pointed that such sequence could be due to selection effects. 

Within the blazar population, only the most powerful objects can be detected at high redshifts 
and they generally belong to the FSRQ class of AGN. According to the ‘blazar sequence’, the 
SED peaks of blazars shift to lower frequencies as their luminosity increases. Since at high 
redshift, only the most powerful objects are expected to be visible, one expects the peaks in their 
SEDs to shift to lower frequencies. Thus, in high redshift blazars the synchrotron hump peaks in 
the sub-mm region and the IC hump peaks in the ~MeV band. In such sources, as the synchrotron 
peak shifts to the mm region, thermal emission from the accretion disk is clearly visible in the 
optical band (e.g., Sbarrato et al. 2013). Accordingly, observing IC peak in the 7 -ray band ( 7 -ray 
detection is the most common blazar fingerprint) is more and more difficult as redshift increases. 
However, in such high redshift blazars, the high energy peak can be observed in hard X-rays due 
to its shift from ~GeV to ~MeV range. In fact, the Burst Alert Telescope (BAT) onboard the 
Swift satellite (Gehrels et al. 2004) has detected blazars up to redshift significantly higher than the 
Fermi-Large Area Telescope (Fermi-LAT) itself (see e.g., Baumgartner et al. 2013). Moreover, 
these high redshift blazars generally harbor more than a billion solar mass black holes at their 
centers (e.g., Sbarrato et al. 2013). 

In this work, a detailed multi-wavelength study of a high redshift blazar S5 0836-1-71 {z = 
2.172; Stickel & Kuehr 1993) is presented. This source has shown multiple episodes of flaring 
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activities since the beginning of Fermi operation. In partieular, a giant 7 -ray outburst was deteeted 
by Ferm/-LAT at the end of 2011 (Ciprini & Dutka 2011). By analyzing the publicly available 
observational data, an exhaustive investigation is performed to understand this peeuliar 7 -ray flare 
and the assoeiated nuclear emission proeesses. Moreover, an attempt is also made to study the 
physieal properties of S5 0836-1-71 during a low activity state, ineluding the first time observations 
from hard X-ray focusing telescope NuSTAR (Harrison et al. 2013). 

This paper is structured as follows. In Seetion 2, a brief colleetion of the main observational 
properties of S5 0836-1-71 is reported. Seetion 3 is devoted to data reduetion procedures followed in 
this work. The results are presented in Section 4 and they are discussed in Section 5. Conclusions 
are outlined in Section 6 . Throughout the work, a ACDM eosmology with the Hubble eonstant 
Ho = 71 km s“^ Mpe“^, = 0.27, and Ha = 0.73 is adopted. 


2. General Properties of S5 0836+71 

The high redshift blazar S5 0836+71 is a radio-bright souree (F 5 ghz = 1-68 Jy; Linford et al. 
2012) and exhibits a flat radio speetrum (a = —0.33; Kuehr et al. 1981). Detection of one-sided 
kilo-parsee seale radio jet is also reported (Hummel et al. 1992). Using the observations from the 
Very Long Baseline Interferometry (VLBI), Otterbein et al. (1998) have revealed the emergenee 
of a new superluminal knot elosely assoeiated with the bright 7 -ray state. From the observations 
at 1.6 and 5 GHz under the VLBI Spaee Observatory Program (VSOP; Lobanov et al. 1998), the 
inner strueture of the jet is found to be of helical shape which is further eonfirmed in a reeent 
study by Perueho et al. (2012). Asada et al. (2010) have proposed a spine-sheath structure of the 
jet using the Very Long Baseline Array observations. 

The optical-UV speetrum of S5 0836+71 is dominated by thermal emission from the 
aeeretion disk (e.g., Ghisellini et al. 2010b). A prominent optieal outburst is also reported by 
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von Linde et al. (1993), wherein the souree was found to be brightened by > 0.5 mag. By 
reprodueing the optieal-UV spectrum with a standard accretion disk model, Ghisellini et al. 
(2010b) have constrained the accretion disk luminosity (Ldisk) and black hole mass (Mbh) of S5 
0836+71 as 2.25 x lO'^^ ergs“^ and 3 x 10® Mq respectively. 

Recent studies focusing on high energy X-ray properties along with the multi-wavelength 
modeling confirmed the blazar nature of S5 0836+71 (e.g., Foschini et al. 2006; Sambruna et al. 
2007; Gianni et al. 2011). It is bright in hard X-rays and included in the 70 months Swift-BAI 
catalog (Baumgartner et al. 2013). It was detected by the Energetic Gamma Ray Experiment 
Telescope onboard the Compton Gamma Ray Observatory (Thompson et al. 1993) and found 
to be a variable emitter of > 100 MeV 7 -rays. Since the launch of Fermi, it is regularly 
monitored by EAT and included in the recently released third Fermz-EAT catalog (3EGE; 

The Eermi-EAT Collaboration 2015). Its 7 -ray spectrum is steep, as expected from a high redshift 
blazar folllowing the ‘blazar sequence’, and is modeled by a logParabola model in the 3EGE 
catalog. Moreover, multiple episodes of 7 -ray flaring from S5 0836+71 are also reported (e.g., 
Ciprini 2011), particularly during late 2011 (Ciprini & Dutka 2011). 

A detailed long term multi-wavelength study of S5 0836+71 has recently been done by 
Akyuz et al. (2013) with a major focus on its variability characteristics. They found the source 
to show flux variations across the electromagnetic spectrum with the strongest variations seen in 
the 7 -ray band. They have not found any significant correlation between 7 -ray and radio fluxes 
and concluded that the observed radio emission can originate from components other than that 
produce 7 -rays. A significant curvature in the flaring state 7 -ray spectrum is also reported by 
them. However, in another study, Jorstad et al. (2013) have proposed the emergence of a new 
radio knot that coincided with the 7 -ray flaring period in 2011 April and they argue that the 7 -ray 
emission region is located at a distance of ~35 pc from the central black hole. 
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3. Multi-wavelength observations 
3.1. Fermi-Large Area Telescope Observations 

The Fermi-LAI data used in this work were eolleeted over the first ~75 months of 
Fermi operation (2008 August 5 to 2014 November 15 or MJD 54,683—56,976). The standard 
data analysis proeedures as mentioned in the Fermi—LAT doeumentation* are adopted. In 
the energy range of 0.1—300 GeV, events belonging to the SOURCE class are used. Good 
time intervals (GTI) are generated by applying a filter expression of “DATA_QUAL>0 && 
LAT_CONFIG==l” and to avoid contamination from the Earth limb 7 -rays, a cut of 100° is also 
applied on the zenith angle. 

Throughout the analysis, the unbinned likelihood method included in the pylikelihood 
library of Science Tools (v9r33p0) and post-launch instrument response functions 
P7REP_SOURCE_V15 are used. Computation of the significance of the 7 -ray signal is performed 
by means of the maximum likelihood (ME) test statistic TS = 2A log(£) where C, represents the 
likelihood function, between models with and without a point source at the position of the source 
of interest. All the sources lying within 10° region of interest (ROI), centered at the position 
of S5 0836-1-71 and defined in the 3EGE catalog, are included for the analysis. The model file 
also includes recently released galactic diffuse emission component gll_iem_v0 5_revl. fit 
and an isotropic component iso_source_v05_revl. txt, as background models^. All the 
parameters except the scaling factor of the sources within the ROI are allowed to vary during 
the likelihood fitting. The normalization parameters of the background models are also left 
free to vary. In addition to that, all the sources lying between 10° to 15° are also considered 
in the analysis. Their spectral parameters are kept fixed to the 3EGE catalog value. A first run 

'http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 

^http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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of the ML analysis is performed over the period of interest and all the sourees with TS < 25 
are removed from the model. This updated model is then used for further analysis. Though S5 
0836+71 is modeled by a logParabola model in the 3FGL eatalog, to generate light eurves, a 
power law (PL) model is used as the PL indiees obtained from this model show smaller statistical 
uncertainties when compared to those obtained from complex model fits. Moreover, since the aim 
is to probe the shortest timescales (hence lower photon statistics), adopting a simple PL model 
is appropriate. For the temporal and spectral studies, the source is considered to be detected if 
TS > 9 which corresponds to ~ 3cr detection (Mattox et al. 1996). Bins with TS < 9 are not 
considered. Primarily governed by uncertainty in the effective area, the measured fluxes have 
energy dependent systematic uncertainties of around 10% below 100 MeV, decreasing linearly in 
log(E) to 5% in the range between 316 MeV and 10 GeV and increasing linearly in log(E) up to 
15% at 1 TeV^. All errors associated with the LAT data analysis are the la statistical uncertainties. 


3.2. NuSTAR Observations 

S5 0836+71 was observed twice with NuSTAR (PI: Fiona Harrison) each in 2013 December 
and 2014 January for a total elapsed time of ~50 and ~70 ksec respectively^. The data in 3—79 
keV are cleaned and filtered for background events using the NuSTAR Data Analysis Software 
(NUSTARDAS) version 1.4.1. Calibration files from NUSTAR CALDB, upadated on 2014 
November 14, are used for instrument responses. All light curves and spectra are extracted for the 
two focal plane modules (FPMA and FPMB) using the tool nuproducts. The source spectra 
are extracted from a 30" circular region centered on S5 0836+71. A circular region of 70" radius 
is selected on the same chip, free from contaminating sources, to extract background spectra. The 

^http://fermi.gsfc.nasa.gov/ssc/data/analysis/LATxaveats.html 

^Due to low Earth orbit, the net exposure time is roughly 50% of the total observation length. 
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source spectra are binned to have at least 20 counts per bin to perform spectral fitting. Further, 
light curves are generated by summing FPMA and FPMB count rates, subtracting background, 
and using 4 ksec binning. 


3.3. Swift Observations 

The Swift satellite has observed S5 0836+71 multiple times with all the three instruments: 
the Burst Alert Telescope (BAT; Barthelmy et al. 2005, 15—150 keV), the X-ray Telescope (XRT; 
Burrows et al. 2005, 0.3—10 keV) and the Ultraviolet Optical Telescope (UVOT; Roming et al. 
2005) which can observe in six filters, namely, V, B, U, UVWl, UVM2 and UVW2. 

The XRT data are first processed with the XRTDAS software package (v.3.0.0) available 
within the HEASOFT package (6.16). Following standard procedures (xrtpipeline 
V. 0.13.0), event files are cleaned and calibrated with the calibration database updated on 2014 
November 7. Standard grade selections of 0—12 in the photon counting mode are used. Energy 
spectra are extracted from the summed event files. Source region is selected as 55" circle while 
background region is chosen as an annular ring with inner and outer radii of 110" and 210" 
respectively, centered at the position of S5 0836+71. Moreover, when the source counts rate 
exceeds 0.5 counts s“^, to avoid pile up effect, the source region is chosen as an annulus with 
inner and outer radii of 5" and 65" respectively, while the background spectra are extracted from 
an annular region of inner and outer radii 130" and 230" respectively (see e.g., Stroh & Ealcone 
2013). The tool ximage is used to combine exposure maps and the ancillary response files are 
generated using the task xrtmkarf. Source spectra are binned to have at least 20 counts per 
bin, using the task grppha and the spectral fitting is performed with XSPEC (Amaud 1996). 
An absorbed power law (A^h = 2.83 x 10^° cm“^; Kalberla et al. 2005) is used for fitting and the 
uncertainties are calculated at 90% confidence level. 
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Swift-\JYOT observations are integrated using uvot imsum and the magnitudes are extraeted 
using the task uvot source. Souree region is seleeted as a eirele of 5" radius eentered at the 
souree position, while the background is chosen from a nearby source-free circular region of 1 ' 
radius. Observed magnitudes are corrected for galactic extinction following Schlafly & Finkbeiner 
(2011) and converted to flux units using the calibrations of Breeveld et al. (2011). These fluxes 
are also corrected for possible absorption by intervening Lyman-a absorption systems, following 
Ghisellini et al. (2010b). 


4. Results 

4.1. Multi-band Temporal Variability 

The long term multi-band light curves of S5 0836-1-71, covering the period 2008 August 5 
to 2014 November 15, are presented in Figure 1. In this plot, Fermi-hNI data points are weekly 
binned and Swift observations correspond to one point per observation ID. The flaring activities 
are seen both in X-ray and 7 -ray bands. In comparison to high energy bands, available optical-UV 
observations show little or no flux variations. This suggests that even during the flaring periods, 
the optical-UV spectrum is accretion disk dominated, which is, in general, not expected to vary 
drastically. To substantiate this, fractional rms variability amplitude (F^ar; Vaughan et al. 2003) is 
calculated and the results are given in Table 1. This parameter is found to be largest for 7 -rays 
and decreases with frequency, a trend generally found among blazars (e.g., Zhang et al. 2005; 
Vercellone et al. 2010). 

It can be seen in Figure 1 that the largest 7 -ray flare, during the period MJD 55,860—55,930, 
is not associated with that in X-rays. The enhancement in X-ray flux is relatively moderate at 
the time of the 7 -ray flare, however, when X-ray flux peaks (MJD 56,790—56,865) there is mild 
increase in the 7 -rays. To study this peculiar behavior, a separate 7 -ray flaring period (Fq; MJD 
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55,860—55,930) and X-ray flaring period (Fx; MJD 56,790—56,865) are seleeted for further 
analysis. For comparison, a low activity state (Q; MJD 56,000—56,657) is also chosen. These 
periods are shown with dashed lines in Figure 1 . The y-ray flaring period selected in this work is 
same as that taken by Akyuz et al. (2013). 

The good photon statistics during the period Fq allows to generate the one day binned y-ray 
light curve, and is shown in the top panel of Figure 2. In order to determine the shortest flux 
variability timescale and the highest 7 -ray flux, the shorter duration flares are further selected from 
the one day binned 7 -ray light curve to generate finer time binned light curves. They are annotated 
by PI, P2, P3, and P4 in Figure 2. As can be seen, the largest flux enhancement is observed during 
the period PI. In the energy range of 0.1—300 GeV, the highest three hours binned 7 -ray flux is 
(5.22 ± 1.10) X 10“® ph cm“^ measured on MJD 55,866, and the associated photon index is 
having a value of 2.62 ± 0.27. This corresponds to an apparent 7 -ray luminosity of (1.62 ± 0.44) 
X 10^° erg s“^. Within error, this luminosity is equal to the highest ever measured from a blazar 
(2.1 ± 0.2 X 10^° erg s“^ for 3C 454.3; Abdo et al. 2011). Further, to have a better estimate of the 
fastest rising and decaying timescales, two rapid flares are selected for time profile fitting. They 
are quoted as FI and F2 in the period PI in Figure 2 and are fitted using a rise and decay function 
given by (see e.g., Abdo et al. 2010a) 



( 1 ) 


where decides the quiescent level and Fq the amplitude of the flares, fp measures the 
approximate time of the peak which is fixed at the time when the highest flux is observed, and 
Tr and Td measure the rise and decay times. The results of the fitting are shown in Figure 3 and 


the obtained parameters are provided in Table 2. The fastest flux rise and decay times are found 


to be ~3 hours and ~2.5 hours respectively. This is the first report of hour scale 7 -ray variability 
detected from this source. Moreover, such short timescale flux variability was also probably not 
known in any high redshift blazar beyond redshift 2 , prior to this work. 
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The maximum X-ray flux from S5 0836-1-71 is measured during the period Fx- In the energy 
range of 0.3—10 keV and eorreeted for galactie absorption, it is having a value of 4.78lQ3g ^ 
10“^^ erg cm“^ s“^, deteeted on MJD 56,836. The assoeiated photon index is hard and has a 
value of 1 . 47 ^ 0 Accordingly, the isotropic X-ray luminosity measured is Lx = 9.43 x 10“^^ 
erg s“^. 

The hard X-ray light curves, using the NuSTAR observations, are generated with 4 ksec 
binning and shown in Figure 4. The flux variations can be seen in both the light curves with the 
variability probability (Abdo et al. 2010a) of 98.33% and 99.86% for 2013 December and 2014 
January observations respectively. Moreover, there is significant change in the flux between the 
two epochs. Earlier, the hard X-ray flux of this source was found to vary on time scales of months 
using the Swift-BAI data (Sambruna et al. 2007). Better sensitivity of NuSTAR has now enabled 
the detection of a shorter time scale of variability from S5 0836-1-71. 


4.2. Spectral Analysis 

The search for the presence of curvature in the 7 -ray spectra is done in all the three 
activity states considered in this work. In particular, two spectral models are used: power law 
(dN/dE (X 77%), where F^.^ is the photon index and logParabola ( dN/dE oc [E /^ 
where Eo is an arbitrary reference energy fixed at 300 MeV, a is the photon index at E^ and fd is 
the curvature index which defines the curvature around the peak). The test statistic of the curvature 
T 5'curve = 2(log E(LogParabola) — log £(power-law)), is calculated to test for the presence of 
curvature. Following Nolan et al. (2012), a threshold TScuwe > 16 is set to test the existence of a 
significant curvature. The resultant fitting parameters are given in Table 3. Significant curvature is 
noticed only during the 7 -ray flare with TS'curve = 51. Though at low significance, there is a hint 
for the presence of curvature in the X-ray flaring period also (TS'curve ~ 9), and during the low 
activity state the power law model represents the spectrum better. 
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To study the overall 7 -ray spectral behavior, the photon index is plotted against the flux in 
Figure 5. In this figure, the variation of photon indices with fluxes is shown both for the weekly 
binned 7 -ray light curve of the entire six years (Figure 1) and the daily binned 7 -ray light curve 
covering the period of GeV outburst (Figure 2). Visual inspection of the daily binned data does not 
reveal any correlation, however, for weekly binned plot the source seems to show a ‘softer when 
brighter’ trend, up to a flux level of ~1.5 x 10“^ ph cm“^ s“^. Above this flux value, there is a 
hint for a ‘harder when brighter’ behavior. To statistically confirm these findings, a Monte Carlo 
simulation is performed that takes into account of the dispersion in the flux and photon index 
measurements. In particular, for each observed randomly drawn pair of flux and photon index 
value, the data is re-sampled by extracting it from a normal distribution centered on the observed 
value and having standard deviation equal to the la error estimate. For daily binned data set, the 
correlation coefficient (p) is 0.02 ± 0.13 with a 95% confidence limit -0.23 ^ p ^ 0.28, clearly 
suggesting the absence of any correlation between fluxes and photon indices. For lower fluxes of 
weekly binned data (i.e., < 1.5 x 10“^ ph cm“^ s“^), p is found to be 0.18 ± 0.12 with 95% 

confidence limit of -0.07 ^ p ^ 0.41. For higher flux level (F.y > 1.5 x 10“^ ph cm“^ s“^), p is 
found to be 0.27 ± 0.11 with 95% confidence limit of 0.02 ^ p ^ 0.48. Thus, based on the Monte 
Carlo analysis, claim for the presence of correlation between fluxes and photon indices in both the 
weekly binned and daily binned data set, as shown in Figure 5, cannot be made. 

Since Swift has monitored S5 0836-1-71 simultaneously with NuSTAR in 2013 December 
and 2014 January (shown by black downward arrows in Figure 1), joint spectral fit to the 
XRT and NuSTAR data is also attempted. The difference in flux calibration between XRT and 
NuSTAR spectra is allowed to vary by including an intercalibration constant (CONST in XSPEC). 
This constant is fixed to 1 for two NuSTAR spectra, since the calibration difference between FPMA 
and FPMB are on the order of 1% (see e.g., Marinucci et al. 2014), and left free to vary for the 
XRT. A first round of fitting is performed with a power law model plus fixed galactic absorption 
which clearly does not provide acceptable fitting parameters (see Table 4; the data/model ratio is 
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shown in Figure 6). However, adding intrinsic absorption (zwabs in XSPEC) or a broken power 
law shape for the continuum, improves the fit (F-test probability of null hypothesis < 10“^°). For 
the case of intrinsic absorption, the level of absorption is higher but with larger uncertainty, and 
the obtained value is similar to that reported by Sambruna et al. (2007). Testing the broken 
power law model on joint XKT-NuSTAR spectra, while keeping the Galactic absorption fixed, also 
results in a satisfactory fit. Moreover, following Akyuz et al. (2013), fitting the data with a power 
law plus two layers of galactic absorption (keeping one absorption fixed to the Galactic value and 
other left to vary) also gives a good fit with statistical parameters similar to that obtained with the 
use of broken power law model (see Table 4). In this case, the additional absorption is higher than 
the Galactic absorption by a factor of ~3. Comparison of the fitting obtained with broken power 
law model with that obtained from power law plus two layers of absorption, are inconclusive. 
The F-test probability of null hypothesis, that the data is better described by the power law plus 
two layers of absorption, is found to be 0.014 and 0.001 for 2013 December and 2014 January 
observations respectively. In Figure 7, the results of fitting the power law with two layers of 
absorption for two simultaneous XKT-NuSTAR observations are shown. 


4.3. Spectral Energy Distributions 

4.3.1. Model Setup 

To reproduce the broadband SEDs of S5 0836-1-71, a simple one zone leptonic emission model 
is developed, by following the prescriptions of Ghisellini & Tavecchio (2009) and Dermer et al. 
(2009) (see also Einke et al. 2008), and is briefly described here. The emission region is assumed 
to be located at a distance of i?diss from the central black hole, and filled with electrons having 
smooth broken power law energy distribution 

(7U~" 

in'HhY + ii'h'^y 


N'ii) = K 


( 2 ) 
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where p and q are the partiele indices before and after the break energy ( 7 (j) respectively 
(primed quantities are measured in the comoving frame). The emission region is assumed to be 
spherical and its size is constrained by considering it to cover the entire jet cross-section with jet 
semi opening angle of 0.1 radians. A standard accretion disk is assumed (Shakura & Sunyaev 
1973) and it produces a multi-temperature blackbody spectrum (Frank et al. 2002). Above and 
below the accretion disk, the presence of X-ray corona is also considered which reprocesses 
30% of the accretion disk luminosity. Its spectrum is considered to be a cut-off power 
law: Lcor(e) oc exp(—e/cc). The BLR is assumed to be a spherical shell located at a 
distance of Rblr = cm, where is the accretion disk luminosity in units of 

10^^ ergs“^ It reprocesses 10% of the accretion disk luminosity. The SED of the BLR is 
approximated as an isotropic blackbody peaking at the rest-frame frequency of the Lyman-a line 
(Tavecchio & Ghisellini 2008). The dusty torus is assumed to be a spherical shell located at a 
distance Rtoms = cm, reprocessing 50% of the accretion disk radiation in the infrared. 

The spectrum of the torus is assumed as a blackbody with temperature Ttorus = epeak’^eC^/S.OSfc, 
where Cpeak is the dimensionless peak photon energy and k is the Boltzmann constant. The 
relative contributions of these emissions with respect to the distance from the central black hole 
are calculated following the methodology of Ghisellini & Tavecchio (2009). In the presence 
of uniform but tangled magnetic field, electrons radiate via synchrotron and SSC mechanisms 
(Finke et al. 2008). The external Compton emissions are also calculated with seed photons coming 
from the accretion disk, the BLR, and the dusty torus (Ghisellini & Tavecchio 2009; Dermer et al. 
2009; Dermer & Menon 2009). The kinetic power of the jet is computed by assuming protons to 
be cold, thus contributing only to the inertia of the jet, and having equal number density to that of 
the relativistic electrons (e.g., Celotti & Ghisellini 2008). 
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4.3.2. SED Modeling Results 

The SED of S5 0836+71 is generated during three different aetivity states. The fluxes are 
averaged for the time intervals shown in Figure 1 and the derived values are given in Table 5. The 
broadband SEDs are generated and modeled using the guidelines presented in Seetion 4.3.1. The 
results are shown in Figure 8 and the assoeiated parameters are presented in Table 6. Variations 
of the energy densities, ealeulated in the eomoving frame, are also shown in the bottom panels of 
Figure 8. 

Two crueial parameters for SED modeling of FSRQs, the aceretion disk luminosity and 
the blaek hole mass, ean be preeisely eonstrained by fitting a standard aeeretion disk model 
(Shakura & Sunyaev 1973) to the optieal-UV speetrum, provided the big blue bump is visible 
(e.g., Ghisellini et al. 2010b). Though this feature is elearly visible in all the aetivity states of S5 
0836+71, the aeeretion disk model is used to reproduee the low aetivity optieal-UV speetrum, as 
during the flaring period there eould be a possibility of eontamination from non-thermal jetted 
emission. The obtained disk luminosity and blaek hole mass are 2.25 x lO'^^ erg s“^ and 3 x 10® 
Mq respeetively, whieh are in perfeet agreement with that obtained by Ghisellini et al. (2010b). 
Aeeordingly, the aeeretion disk luminosity is ~60% of Eddington luminosity. 

Although the observed y-ray variability timeseales are a fraetion of day, longer integration 
times are required to produee good quality EAT SEDs for speetral modeling. The fluxes at other 
frequeneies are also averaged during the y-ray integration period. Therefore, parameters obtained 
from the generation and modeling of the SEDs should be eonsidered as representative of the 
average eharaeteristies of S5 0836+71 during the various aetivity states. 
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5. Discussion 

In this work, major emphasis is given on the 7 -ray temporal properties of S5 0836-1-71 during 
its GeV outburst in 2011 and on the physieal properties of the souree via SED modeling. The 
results of the monitoring from the sensitive hard X-ray teleseope NuSTAR are also presented for 
the first time. 

The long term multi-wavelength study of S5 0836-1-71 is done in great detail by Akyuz et al. 
(2013). They report a weak optieal- 7 -ray flux eorrelations. In this work also, though data is too 
sparse, visual inspeetion suggests the absenee of any analogy between 7 -ray and optieal fluxes. 
This eould be primarily due to the faet that optieal-UV speetrum of S5 0836-1-71 is dominated 
by the aeeretion disk radiation, whieh is, in general, unlikely to vary on short timeseales. An 
interesting feature of the multi-frequeney light eurve is the variation of X-rays as eompared 
to 7 -rays. Though the fluxes have enhaneed during different flaring periods, the amplitude of 
variability does not follow the same pattern (Figure 1). At the time of GeV outburst. X-ray 
flux variations are relatively moderate, and when the X-ray flux is highest, there is relatively 
mild enhaneement in 7 -rays. This is explained by loeating the emission region at different jet 
environments where various radiative energy densities play major role in deseribing the data, as 
diseussed later in this seetion. 

The weekly binned 7 -ray light eurve, presented in the top panel of Figure 1 , show multiple 
episodes of flaring aetivities with the brightest flare oeeurred during 2011. Apart from that, the 
souree is sporadieally deteeted throughout the Fermi operation. The good photon statisties during 
the highest 7 -ray aetivity period makes possible to study this event in detail. A daily binned 7 -ray 
light eurve is generated whieh shows the presenee of sub-flares within the time period of enhaneed 
emission (Figure 2). From the one day binned light eurve, further periods of high aetivity are 
seleeted (PI, P2, P3, and P4) and elose-up light eurves are generated with finer time bins. As ean 
be seen, the flares have eomplex strueture and asymmetrie profiles and in few eases signifieant 
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variations are visible. Results of the flare profile fitting indieate that the flare FI seems to have a 
symmetric profile and the flare F2 possess asymmetric shape with fast rise and slow decay. Such a 
fast rise of the flare could be associated with the faster acceleration, probably at shock front, and 
the slow decay can be attributed to the weakening of the shock. 

The location of the emission region, obtained from the SED modeling, is found to be inside 
the BLR where the energy density provided by the BLR clouds should dominate the energy 
densities of other photon fields in the jet rest frame, and hence the primary mechanism for the 
production of 7 -rays would be IC upscattering of the BLR photons (observed energies eo 
10 eV). Accordingly, the cooling timescale for the electrons emitting 7 -ray radiation with the 
energies e.^ ~ 1 GeV, measured in the observer’s frame, would be 



(3) 


i.e., ~9 minutes for the comoving BLR photon energy density — 10 erg cm“^. The 
calculated cooling time is significantly shorter than the observed shortest flux decay time of the 
flares, implying that the observed flux decrease is governed by processes other than radiative 
energy losses, probably a combination of different factors such as particle acceleration or jet 
dynamics (e.g., Bottcher & Principe 2009; Kushwaha et al. 2014). On the other hand, geometry 
and presence of sub-structures in the emitting region could also led to such observations 
(Tanihata et al. 2001). 

At the peak of the 7 -ray flare, the highest measured isotropic 7 -ray luminosity is ~1.64 x 
Iflso gj-g which corresponds to the total power emitted in the 7 -ray band, in the proper frame 

of the jet, — L^/2V‘^ ~ 2.27 x lO'^^ ergs“\ assuming the bulk Lorentz factor F = 19 
obtained from the modeling of the 7 -ray flaring SED. This is a substantial fraction of the kinetic 
jet power (~15 %, Pj,kin = 1-48 x 10^® ergs“k), implying the jet becomes radiatively efficient 
and the bulk of the radiative energy released in the form of 7 -rays. Lor a comparison, the em 
is found to be about 60% of the Eddington luminosity. Moreover, em is also a good fraction of 
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the entire available aeeretion power (Lace — -^disk/'^disk = 2.25 x 10^® erg s assuming radiative 
effieieney r^disk = 10 %). 

The results on energetics suggest that both the total jet power (dominated by protons), and 
Lace exceed the Eddington luminosity. However, there could be few possible caveats. For a given 
value of Ldisk, a relatively higher r^disk will reduce the total accretion power. Presence of more than 
one leptons per proton in the jet will decrease the budget of the proton jet power. Most importantly, 
underestimation of the central black hole mass can also lead to super Eddington accretion and jet 
power. The mass of the central black hole used in this work is 3 x 10® Mq, which is obtained by 
fitting the standard accretion disk model to the optical-UV spectrum. However, if the C IV line 
parameters from Torrealba et al. (2012) are used to calculate the black hole mass of S5 0836+71 
and the virial relations of Shen et al. ( 2011 ) are adopted, the estimated black hole mass is as large 
as ~1.3 X 10^® Mq. Therefore, these uncertainties and/or their combinations could lead to the 
total jet power and the accretion power exceeding the Eddington luminosity. 

Curvature in the 7 -ray spectrum is a characteristic property of powerful FSRQs (see e.g., 
Abdo et al. 2010b; Harris et al. 2014). Interestingly, this feature is seen to be more prominent 
during the 7 -ray flaring activities (see e.g., Abdo et al. 2011; Paliya et al. 2015, for 3C 454.3 and 
3C 279 respectively). During the period of 7 -ray outburst, a significant curvature is noticed in 
the 7 -ray spectrum of S5 0836+71 (Table 3) similar to that reported by Akyuz et al. (2013). A 
possible explanation of such curvature could be due to the attenuation of 7 -rays by photon-photon 
pair production on He ll Eyman recombination lines within the BER (Poutanen & Stern 2010). 
Moreover, Cerruti et al. (2013) have proposed an alternative model in which a log parabolic 
particle energy distribution and Klein-Nishina (KN) effect on IC scattering of BER photons 
reproduce the observed curvature. Here, the spectral curvature is explained on the basis of KN 
mechanism and broken power law electron distribution. Therefore, though the high redshift and 
steep 7 -ray spectrum suggests that S5 0836+71 would be weak in the 7 -ray window of the SED, 
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the giant GeV outburst of this source reveals many peculiar characteristics which are generally 
seen only in low redshift blazars. 

The data to model ratio of Swift observations, shown in Figure 6 along with the NuSTAR data, 
show a clear hardening of the X-ray speetrum below 1 keV, equally well fitted by a power law 
model with absorption in exeess of the Galactic value or by a broken power law eontinuum with 
Galactic absorption. In fact. X-ray observations of high redshift quasars indieate for the presenee 
of significant excess column densities in radio-loud sources (e.g.. Page et al. 2005; Tavecchio et al. 
2007, and referenees therein). Moreover, in several high redshift blazars the hardening of the 
soft X-ray spectrum is explained in terms of absorption by warm plasma (iVn 10^2 cm- 2 ) 
surrounding the quasar, and probably in the form of outflow/wind (e.g., Fabian 1999). Since 
blazars typically do not show any significant intrinsic absorption, as the jet is expected to remove 
the gas from its vieinity, the presence of large gas column along the line of sight is highly 
improbable. One of the alternative explanation could be the underestimation of Galaetie hydrogen 
eolumn density by a faetor of obtained in this work, similar to that noted by Akyuz et al. 
(2013). In the souree frame, the Fe-K emission line is located at around 6.4 keV, and for S5 
0836-1-71 it is expeeted to be around 2 keV in the observer’s frame. From Figure 6 , it can be seen 
that there is slight excess emission at ~2-3 keV, however, due to short exposure of the XRT and 
hence low eounting statisties, a strong claim cannot be made. 

The homogeneous one zone synchrotron plus IC model with the inelusion of thermal 
emission from the accretion disk, sueeessfully reproduces the observed SEDs. However, it does 
not explain the radio data, suggesting that much of the radio emission does not come from the 
eompact emission region, primarily due to synehrotron self absorption proeess. The observed 
radio emission can be due to a superposition of many jet components (Konigl 1981), which is 
also supported by the lack of correlation seen in radio and 7 -ray fluxes (Akyuz et al. 2013). The 
optieal-UV part of the SEDs in all the aetivity states is dominated by the radiation from the 
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accretion disk and this explains the laek of signifieant variability seen in the optieal-UV band. 
Optieal polarization observations during the peak of the 7 -ray flare reveals the deteetion of very 
small polarization varying between 1% to 6 % (Jorstad et al. 2013), strengthening the idea of the 
dominanee of optieal-UV emission from the aeeretion disk. Further, the X-ray to high energy 
7 -ray window of the SEDs is explained via IC seattering of the photons originating outside the jet, 
primarily from the aeeretion disk and the BLR. Though in many FSRQs the origin of X-rays is 
attributed to SSC process (e.g., Cerruti et al. 2013), it is better explained by EC meehanism for the 
ease of high redshift blazars (Ghisellini et al. 2010b). In this work, the X-ray speetrum is mainly 
reproduced by EC-disk and 7 -ray emission via EC-BLR meehanism, as ean be seen in Figure 8 . 

To eompare the physical properties of S5 0836-1-71 in different activity states, all the three 
SEDs are plotted together in Figure 9. Looking at the obtained parameters in Table 6 , the emission 
region is found to be inside the BLR in all the aetivity states, similar to that found in earlier studies 
(Sambruna et al. 2007; Taveeehio et al. 2007; Celotti & Ghisellini 2008; Ghisellini et al. 2010b). 
Comparing the parameters obtained from the modeling of the low aetivity state SED with that 
obtained during the 7 -ray flare suggests the enhancement of bulk Lorentz factor as a major cause 
of the flare. However, it should be noted that these results ean also be reprodueed by keeping the 
bulk Lorentz faetor fixed and varying the total jet power and the equipartition condition. In such a 
situation, the required jet power would be less for the low activity state. Moreover, there are few 
other ehanges, sueh as slight hardening of the speetral shape of the partiele energy distribution, 
loeation of the emission region (and henee its size), and deerease in the magnetie field. The 
Compton dominanee (whieh is the ratio of IC to synehrotron peak luminosities) also inereased by 
about an order of magnitude. Though the synehrotron flux also enhaneed during the 7 -ray flare, 
its eontribution to the observed optieal-UV speetrum is negligible and thus no major ehange in the 
optical band is noted. 


While the three SED models are elose to equipartition between eleetrons and magnetie fields 


- 22 - 


{Pe rsj Pb), the magnetic power is tiny compared to the proton power, or total jet power. This 
suggests a very low magnetization of the emitting region. Moreover, the radiative efficiency of 
the jet (expressed as Pj-Z-fjet), for Fq, Q, and Fx periods, are obtained as 0.23, 0.03, and 0.20 
respectively. This suggests that during the flares, a good fraction of the total jet power gets 
converted to the radiative power. However, majority of the jet power remains in the form of kinetic 
power and used to produce large scale jets. Looking at the various jet powers in Table 6 , it is clear 
that both Pe and Pb are smaller than P^. The fact that Pg < Pr appears to be a contradiction at 
first sight since the radiation produced by electrons, in principle, cannot have more power than 
the emitting electrons. However, this inference may not be correct because the cooling timescale 
of 7 -ray producing electrons are much shorter than Pbiob/c. Therefore, in order to maintain the 
emission over the entire flare period, the number of high energy electrons should be replenished 
till the flare duration. In other words, the electron population continuously gets energized at the 
expense of bulk energy and transfer the power to radiation. Hence Pj. should be treated as the 
fraction of bulk luminosity transferred to radiation through electrons. On the other hand, since the 
mean energy of the electrons will be closer to 7 min, Pe is mainly decided by the total number of 
electrons at lower energies which do not emit at 7 -rays. 

An interesting phenomenon is the detection of large X-ray flare with relatively moderate 
7 -ray variability amplitude (see period Fx in Figure 1). Since in this work, the X-ray emission 
is mainly explained by EC-disk mechanism, this suggest that a possible cause of the X-ray flare 
could be enhancement in the EC-disk flux. One of the possibility, perhaps, is the emission 
region to be located close to the central black hole. This is because the energy density of the 
direct disk radiation, as seen in the comoving frame, increases with decrease in the dissipation 
distance (e.g., Dermer et al. 2009; Ghisellini & Tavecchio 2009). The variation of comoving 
frame energy densities as a function of the distance of the emission region from the central black 
hole for different activity periods, are shown in the bottom panels of Eigure 8 . It can be seen in 
this plot that during the X-ray flaring period, the emission region is located in a jet environment 
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where the accretion disk energy density dominates over other energy densities. Due to this effect, 
the enhancement of EC-disk radiation will be higher as compared to EC-BLR component and 
accordingly the rise in the X-ray flux will be more. 


6. Summary 

The detailed multi-frequency study of the high redshift blazar S5 0836-1-71 is presented in 
this work. The main findings are as follows 

1. The long term multi-wavelength light curves of S5 0836-1-71 show more than one episodes 
of flaring activities in X-ray and 7 -ray bands, while optical-UV fluxes show little or no 
variations. 

2. The highest three hours binned 7 -ray flux and associated photon index are found to be (5.22 
± 1.10) X 10“® ph cm“^ and 2.62 ± 0.27 respectively, measured on MJD 55866. This 
corresponds to an isotropic 7 -ray luminosity of (1.62 ± 0.44) x 10^° erg s“\ thus making it 
one of the most luminous 7 -ray flare ever observed from the blazar class of AGN. 

3. The observed 7 -ray flux variability of about three hours is the shortest flux rising/decaying 
time ever measured from this source, and probably from any high redshift blazar beyond 
redshift 2 . 

4. A search for possible correlation between 7 -ray flux and photon index, using a Monte Carlo 
simulation approach, does not result in any significant correlation between them. 

5. Reproduction of the observed SEDs with a simple one zone leptonic emission model 
suggests the emission region to be located inside the BER, and the optical-UV spectrum is 
dominated by the accretion disk radiation. 
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6 . Enhancement in the bulk Lorentz faetor is probably a primary faetor of the observed giant 
7 -ray outburst. The high aetivity seen in the X-ray band with less variable 7 -ray eounterpart 
ean be explained by locating the emission region closer to the eentral blaek hole where 
the eomoving frame energy density from the accretion disk dominates over other energy 
densities, resulting in enhaneed EC-disk flux peaking in X-rays. 
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Table 1: Fractional rms variability amplitude parameter (F^ar), ealeulated for the light eurves shown 
in Figure 1 and 4. 


Energy band 

F 

^ var 

7 -ray (0.1-300 GeV) 

0.79 ± 0.02 

X-ray (3-79 keV) 

0.21 ± 0.01 

X-ray (0.3-10 keV) 

0.33 ± 0.01 

UVW2 

0.07 ± 0.02 

UVM2 

0.04 ± 0.05 

UVWl 

0.09 ± 0.02 

u 

0.05 ± 0.02 

B 

0.04 ± 0.03 

V 

0.04 ± 0.08 


Table 2: Parameters obtained from the time profile fitting of two flares. Errors are estimated at la 
level._ 


Name 

Fc 


tp 

Tr 

Td 

X? 

El 

1.11 ±0.27 

8.41 ±2.01 

55866.31 

0.13 ±0.05 

0.10 ±0.04 

0.77 

E2 

1.31 ±0.26 

3.54 ± 1.29 

55872.75 

0.08 ± 0.06 

0.27 ±0.13 

0.21 


Note. — Fluxes Fc and Fp are in 10 ® ph cm ^ s Tr and Td are in days, and tp, which is kept fixed, is in MJD. 
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Table 3: Results of the Model Fitting to the 7 -ray Spectra of S5 0836-1-71, obtained for differ¬ 
ent time periods. Col.[l]: period of observation (MJD); Col.[2]: activity state; Col.[3]: model 
used (PL: power law, LP: logParabola); Col.[4]: integrated 7 -ray flux (0.1—300 GeV), in units of 
10“^ ph cm“^ s“^; Col.[5] and [ 6 ]: spectral parameters (see definitions in the text); Col.[7]: test 
statistic; Col.[ 8 ]: test statistic of the curvature. 


Period 

[1] 

Activity 

[2] 

Model 

[3] 

-fo.l-300 GeV 

[4] 

Fo.l-soo Gev/a 

[5] 

/3 

[6] 

TS 

[7] 

T 5*curve 

[8] 

55,860-55,930 

7 -ray flare (Fq) 

PL 

6.25 ± 0.21 

2.67 ± 0.04 


2621.34 




LP 

5.85 ± 0.22 

2.50 ± 0.05 

0.29 ± 0.05 

2638.81 

51.0 

56,000-56,657 

Low activity (Q) 

PL 

0.51 ± 0.04 

3.09 ±0.10 


242.03 




LP 

0.49 ± 0.04 

3.04 ±0.11 

0.13 ±0.12 

241.75 

1.26 

56,790-56,865 

X-ray flare (Fx) 

PL 

2.22 ±0.15 

2.59 ± 0.06 


582.70 




LP 

2.05 ±0.16 

2.40 ±0.10 

0.18 ±0.07 

586.40 

8.96 
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Table 4: Summary of the Joint XKT+NuSTAR Spectral Fitting. Col.[l]: fitted model (PL: power 
law, BPL: broken power law, GA: galactic absorption, AA: additional galactic absorption, zA: 
absorption in the intrinsic source frame); Col.[2]: additional neutral Hydrogen column density 
(10^° cm~^); Col.[3]: intercalibration constant ; Col.[4]: photon index of PL model or photon 
index before break energy in BPL model; Col.[5]: photon index after break energy in BPL model; 
Col.[6]: break energy (keV); Col.[7]: normalization at 1 keV (10“^ ph cm~^ s“^ keV~^); Col.[8]: 
statistical parameters. 


Model 

^addi. 

CONST 

Tx/ri 

Ta 

-^break 

Norm. 

X^/dof. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 





2013 December 15 




PL+GA 


0.81 ±0.03 

1.58 ±0.01 



3.02 ± 0.09 

692.30/550 

PL+GA+zA 

87.25 ± 10.94 

0.85 ± 0.03 

1.64 ±0.02 



3.48 ±0.12 

573.44/549 

PL+GA+AA 

8.31 ±0.93 

0.88 ± 0.03 

1.66 ±0.02 



3.57 ±0.12 

559.43/549 

BPL+GA 


0.93 ± 0.06 

1.13 ±0.10 

1.67 ±0.03 

2.07 ± 0.37 

2.46 ±0.18 

553.30/548 





2014 January 18 




PL+GA 


0.72 ± 0.02 

1.66 ±0.01 



5.54 ±0.13 

880.74/762 

PL+GA+zA 

106.51 ± 13.84 

0.83 ± 0.03 

1.69 ±0.01 



5.90 ±0.14 

779.32/761 

PL+GA+AA 

9.58 ± 1.11 

0.87 ± 0.03 

1.69 ±0.01 



5.97 ±0.14 

766.55/761 

BPL+GA 


0.98 ± 0.08 

1.23 ±0.10 

1.70 ±0.02 

2.83 ± 0.62 

3.71 ±0.33 

755.79/760 
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Table 5: Summary of the SED analysis. Fermi-hKY analysis results are given in Table 3. 

NuSTAR 


Activity state 

Exp.®^ 

T t) 

i 3-79 keV 

F 3-79 keV'^ 

Norm.'^ 

Stat.'^ 


Q 

29.7 

1.64 ± 0.04 

4 99+0-19 
+•^^-0.21 

3.38 ± 0.28 

425.79/436 





Swift-XRT 




Activity state 

Exp.'^ 

Fq.S-IO keV^ 

-fo. 3-10 keV® 

Norm.'^ 

Stat.*^ 


Eg 

20.65 

1.21 ±0.02 


2.58 ± 0.06 

265.68/313 


Q 

8.06 

1.34 ±0.04 

2 36+0-09 

+■-’'9-0.10 

2.35 ± 0.09 

143.15/126 


Fx 

2.32 

1.26 ±0.06 

4 37+0.28 
^■■’'-0.28 

3.98 ± 0.24 

80.88/56 





Swift-UVOT 




Activity state 

yh 


\J^ 

UVWl*’ 

UVM2’’ 

UVW2*’ 

Fg 

3.73 ±0.16 

4.74 ±0.14 

5.66 ±0.12 

3.57 ± 0.05 

3.31 ±0.05 

3.20 ± 0.03 

Q 

3.66 ± 0.20 

4.72 ±0.18 

5.62 ±0.18 

3.60 ±0.11 

3.47 ± 0.06 

3.16 ±0.05 

Fx 

3.73 ± 0.33 

4.48 ± 0.25 

5.78 ±0.18 

3.68 ±0.13 

3.10 ±0.10 

3.16 ±0.06 


^Net exposure in kiloseconds. 

^Photon index of the power law model in 3—79 keV energy range. 

^Power law flux in 3—79 keV energy range, in units of 10“^^ erg cm“^ s“^. 
^Normalization at 1 keV in 10“^ ph cm“^ keV~^. 

^Statistical parameters: x^/dof. 

®Photon index of the absorbed power law model in 0.3—10 keV energy band. 
^Unabsorbed flux in units of 10~^^ erg cm“^ s“^, in 0.3 — 10 keV energy range. 
^Average flux in Swift V, B, U, Wl, M2, and W2 bands, in units of 10“^^ erg cm“^ s“^. 
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Table 6: Summary of the parameters used/derived from the modeling of the SEDs in Figure 8. 
Viewing angle is taken as 3° and the eharaeteristie temperature of the torus as 900 K. For a disk 
luminosity of 2.25 x 10^^ erg s“^ and black hole mass of 3 x 10® Mq, the size of the BFR is ~0.5 
parsec (1700 i?sch)- 




Activity state 


Parameter 

Fg 

Q 

Fx 

Slope of particle spectral index before break energy (p) 

1.6 

1.7 

1.6 

Slope of particle spectral index after break energy (q) 

4.0 

4.2 

3.8 

Magnetic field in Gauss (B) 

1.7 

1.8 

2.5 

Equipartition factor* (r/equi.) 

0.89 

2.74 

1.43 

Bulk Lorentz factor (F) 

19 

10 

15 

Break Lorentz factor 

84 

69 

73 

Maximum Lorentz factor (7max) 

5e4 

5e4 

2e4 

Dissipation distance in parsec (Rsch) 

0.14 (500) 

0.11 (400) 

0.07 (250) 

Jet power in electrons in log scale (Pe) 

45.83 

45.62 

45.57 

Jet power in magnetic field in log scale (Pb) 

45.88 

45.18 

45.41 

Radiative jet power in log scale (P,-) 

47.84 

46.58 

47.23 

Jet power in protons in log scale (Pp) 

48.17 

48.06 

47.93 


Note. — *Equipartition factor is the ratio of the particle to magnetic energy density (pequi. = F'/C/g). 
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Fig. 1.— Multi-frequency light curves of S5 0836-1-71 since the launch of Fermi satellite. Fermi- 
LAT and Swift-XKT data points are in units of 10“^ ph cm”^ s“i and counts s“^ respectively. UV 
and optical fluxes have units of 10“^^ erg cm”^ s“i . Periods corresponding to high 7 -ray and 
X-ray states, and a low activity phase are quoted as Fq, Fx, and Q respectively. Black downward 
arrows represent the time of NuSTAR observations. See text for details. 
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Fig. 2.— Gamma-ray light curves of S5 0836-1-71 covering the period of GeV outburst. Fluxes are 
in units of 10“® ph cm“^ s“^. Good photon statistics during this period allows to go for various 
time binnings. Top panel represents the daily binned light eurve, from whieh, intervals of higher 
aetivities (marked as PI, P2, P3, and P4) are selected to go for twelve hours, six hours, and three 
hours binning. In the three hour binned light eurve of PI period (bottom left eorner), two flares FI 
and F2 are seleeted for prohle htting. 
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Fig. 3.— Top: Three hours binned 7 -ray flares of S5 0836-1-71 seleeted for the time profile fitting, 
during the period of GeV outburst. Fluxes have same units as in Figure 2. Blue solid line denotes 
the best-fit temporal profile assuming an exponential rise and fall. Bottom: The residual of the 
fitting. 
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Fig. 4.— Baekground subtraeted 3—79 keV light eurves of the two NuSTAR observations of S5 
0836-1-71. For this plot, the FPMA and FPMB count rates are summed, and 4 ksec binning is 
applied. 
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Fig. 5.— Left: Weekly seatter plot of the 7 -ray flux vs. photon index. Right: Same as left but for 
daily binned data eovering the period of high 7 -ray aetivity. 
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Fig. 6 .— Data to model ratio of the fit of joint XRT and NuSTAR speetra with a power law and 
fixed galaetie absorption. The defieit of soft photons below 1 keV elearly indieates a greater degree 
of absorption or an intrinsic hardening of the spectrum. 
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Fig. 7.— Top: Joint XRT (0.3—10 keV) and NuSTAR (3—78 keV) spectrum, fit with power law 
plus two layers of absorption. Bottom: data to model ratio for the obtained fit. Two NuSTAR spectra 
are grouped in XSPEC and rebinned for visual clarity. 
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Fig. 8.— Top: Spectral energy distributions of S5 0836+71 during different activity states. Si¬ 
multaneous Swift and Fermz-LAT data are shown with red squares and light blue circles represent 
archival observations. During the low activity period, the NuSTAR data is also available which is 
shown by green squares. Pink thin solid and green long dashed lines correspond to synchrotron 
and SSC respectively. EC-disk and EC-BLR mechanisms are shown by grey dashed and orange 
dash-dash-dot lines respectively. Blue thick solid line is the sum of all the radiative contributions. 
Bottom: Variations of the energy densities measured in the comoving frame, as a function of the 
distance from the black hole, in units of i?sch- Vertical line denotes the location of the emission 


region. 
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Fig. 9.— The SEDs of S5 0836+71 in different aetivity states, same as in Figure 8, but plotted 


together for comparison. 







